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Inflationary cosmology provides a natural mechanism for the generation of primordial perturba- 
tions which seed the formation of observed cosmic structure and lead to specific signals of anisotropy 
in the cosmic microwave background radiation. In order to test the broad inflationary paradigm as 
well as particular models against precision observations, it is crucial to be able to make accurate 
predictions for the power spectrum of both scalar and tensor fluctuations. We present detailed 
calculations of these quantities utilizing direct numerical approaches as well as error-controlled uni- 
form approximations, comparing with the (uncontrolled) traditional slow-roll approach. A simple 
extension of the leading-order uniform approximation yields results for the power spectra ampli- 
tudes, the spectral indices, and the running of spectral indices, with accuracy of the order of 0.1% 
- approximately the same level at which the transfer functions are known. Several representative 
examples are used to demonstrate these results. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Cosmological inflation Q is a central component of 
the present theoretical picture of cosmology. Inflation 
not only directly addresses and solves fundamental weak- 
nesses of the older Big Bang picture - the flatness and 
horizon problems - it also provides an elegant mechanism 
for the creation of primordial fluctuations Q, essential 
for explaining the observed structure of the present-day 
Universe. As the Universe inflates, microscopic quantum 
scales are stretched to macroscopic cosmological scales. 
Quantum vacuum fluctuations provide the initial seeds 
that are amplified by gravitational instability, leading to 
the formation of structure in the Universe. 

The primordial fluctuations associated with inflation 
are, primarily, of a very simple type. In "standard" in- 
flationary models, they arise from the fluctuations of an 
effectively free scalar field and are hence Gaussian ran- 
dom fields, completely characterized by two-point statis- 
tics, such as the power spectrum. The basic task, then, 
is to determine the scalar and tensor perturbations (vec- 
tor perturbations being naturally suppressed) in terms of 
the associated power spectra. Power spectra from stan- 
dard inflation models are conveniently parameterized by 
a spectral index and its (weak) variation with scale - the 
"running" of the spectral index. In addition, inflation 
couples the generation of scalar and tensor fluctuations; 
if both are independently measured, "compatibility" re- 
lations characteristic of inflation can be put to observa- 
tional test. It should be noted that complicated infla- 
tionary models can be constructed to produce baroque 
"designer" power spectra. Present observations do not 
require building such models. 

Observational constraints on the primordial power 
spectra arise from measurements of temperature 
anisotropies in the cosmic microwave background radia- 



tion (CMBR), sensitive to both scalar and tensor modes, 
and from measurements of the density power spectrum 
in surveys of the large-scale mass distribution in the Uni- 
verse. Ground and satellite-based observations of the 
CMBR anisotropy have yielded results very consistent 
with the essential adiabatic and Gaussian nature of in- 
flationary perturbations, with a value of the scalar spec- 
tral index very close to unity |3| . The tensor fluctuations 
are expected to be much lower in amplitude and have 
not yet been observed. Scalar perturbations seed struc- 
ture formation and hence can be measured by observing 
the large-scale distribution of galaxies and neutral hydro- 
gen, as in ongoing redshift surveys 3 and Lya observa- 
tions [U- This second set of independent measurements 
provides information on smaller scales than the CMBR, 
yet there is an overlap region where both measurements 
have been shown to be consistent (Fig. [I}. The observed 
power spectrum is P bs(k) — P(k)T 2 (k) where P{k) is 
the primordial power spectrum and T(k) is the transfer 
function for radiation or matter as appropriate. 

As measurements continue to improve, tests of the in- 
flationary paradigm and specific inflationary models will 
become more stringent, especially if accurate measure- 
ments of the running of the spectral index are performed 
and tensor perturbations are observed. In addition, an 
inverse analysis of the observational data may then be at- 
tempted, in an effort to directly measure the inflationary 
"equation of state" 0. In keeping with the remarkable 
improvement of observational capabilities, the quality of 
theoretical predictions has necessarily to be addressed. 
The central theme of this paper is to quantify and im- 
prove the quality of theoretical predictions. 

Scalar and tensor power spectra for particular infla- 
tion models can be calculated either by direct numerical 
methods or by employing analytic or semi-analytic ap- 
proximations such as the slow-roll approach 8] or uni- 



2 



Wavelength X [h _1 Mpc] 
10" 1000 100 10 1 




Wavenumber k [h/Mpc] 



FIG. 1: Compilation of observations of the matter power spec- 
trum Pobs(fc) taken from Tegmark et al. || (with permission). 
This particular figure assumes a flat scalar scale-invariant 
model with Q m = 0.28, h = 0.72 and Q, b /Q. m = 0.16, r = 0.17, 
and a bias, 6, = 0.92 for the Sloan Digital Sky Survey (SDSS) 
galaxies. The solid line is the theoretical curve. 



form approximations 0, . For individual models, the 
numerical approach may well be preferable, but powerful 
approximations have their own advantages. They provide 
intuition and understanding applicable to entire classes of 
models. Provided tight error controls can be met, they 
are much faster than a mode-by-mode numerical inte- 
gration for obtaining the power spectrum. Also, both 
numerical and approximate strategies can be melded to- 
gether by first obtaining a set of approximate results over 
a wide range of parameters and then spot-checking nu- 
merically. It is our view that both strategies should be 
applied and compared against each other as they have 
differing error modes. Finally, it should be kept in mind 
that the ultimate accuracy with which the primordial 
fluctuations must be computed is limited by how accu- 
rately the radiation and matter transfer functions are 
known, as well as by the errors associated with observa- 
tions. At present, the transfer functions can be computed 
to about 0.1% accuracy 

.11]. 

sufficient for dealing with 
next-generation observations. We will use this value as a 
relative figure of merit when discussing errors below. It 
should be noted that the accuracy requirements for the 
tensor or gravitational wave component are not as strin- 
gent as for the scalar component; the expected signal 
is at large length scales, where cosmic variance becomes 
unavoidable. Finally, measurements of the radiation and 
matter power spectrum at higher values of k face issues 
such as the rapid fall-off of the primary signal, secondary 



point-source contamination, foreground subtraction (for 
the CMBR), galaxy bias, systematic errors, effects from 
baryons and neutrinos 0, and limited accuracy in 
theoretical computation of the present nonlinear power 
spectrum [l4j for the matter distribution. 

In this paper we focus on the perturbation spectrum for 
single- field inflationary models. We implement a direct 
numerical approach as well as the slow-roll and uniform 
approximations, with the aim of establishing control over 
the errors associated with each method and understand- 
ing their associated advantages and disadvantages. The 
main technical advances are a robust and numerically ef- 
ficient strategy for the mode-by-mode integration, simple 
and useful error estimates for the uniform approximation, 
and a simple improvement strategy for power spectra am- 
plitudes for the uniform approximation at leading order. 
We are able to show that the improved leading order uni- 
form approximation leads to very good accuracy for the 
spectral indices and their running, as well as for the am- 
plitudes of the power spectra and the ratio of tensor to 
scalar perturbations. For the most part, the accuracy of 
the results obtained is of the order of 0.1%. 

The method of uniform approximation employed in 
this paper has been presented in Refs. H [[(I. This 
method is a "uniformization" of the well-known Wentzel- 
Kramers-Brillouin (WKB) or Liouville-Green (LG) ap- 
proximation |15| in the presence of transition points. 
The uniform approximation began with the work of 
Langer an d others, and was followed by the notable 
contributions of Olver [rH fl8| which we rely on for our 
analysis below. While this line of investigation rested on 
the analysis of ordinary differential equations, equivalent 
asymptotic results based on an integral representation 
were given by Chester, Friedman, and Ursell The 
uniform approximation has proven very useful in many 
applications, e.g., chemical physics |2fJ, the semiclassi- 
cal limit in quantum mechanics |2l| . and the quantum- 
classical transition in quantum cosmology |22j |. For our 
purposes, the key advantages of the Olver uniform ap- 
proximation are that it does away with WKB-like match- 
ing conditions j2^| (such a procedure fails to indicate the 
error of the approximation) , has controlled error bounds 
over the entire domain of interest, is systematically im- 
provable, and allows analytic simplifications in special 
cases of physical interest. 

The paper is organized as follows. Section ITT1 provides 
the necessary background regarding the calculation of the 
primordial power spectra and spectral indices for single 
field inflation models. In Section llTTl the results from the 
uniform approximation in leading order are given, while 
in Section Hvl the essential equations for the slow-roll ap- 
proximation are summarized. In Section we give a de- 
tailed description of the numerical determination of the 
primordial power spectra and spectral indices by solving 
the mode-equations numerically, in the uniform approx- 
imation, and in the slow-roll approximation. We investi- 
gate three different examples in Section IVTl and conclude 
with a discussion of our results in Section lYlII 
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II. BACKGROUND 

The generation of perturbations during inflation is due 
to the amplification of quantum vacuum fluctuations by 
the dynamics of the background spacetime. In order to 
calculate the perturbation spectra for these quantum vac- 
uum fluctuations three steps are necessary: the dynam- 
ics of the background spacetime must be determined, the 
mode equations for the scalar and tensor perturbations 
must be solved, and finally, the power spectra themselves 
must be calculated as functions of wave-number k. We 
will briefly outline these three steps in the following; for 
detailed derivations the reader is refered to the litera- 
ture |5l2ll2ll2^|. 



B. Mode Equations, Power Spectra, Spectral 
Indices, and their Running 

The modern understanding of fluctuations generated 
by inflation is based on the gauge-invariant treatment 
of linearized fluctuations in the metric and field quanti- 
ties |3 E3 HE HH ■ A particularly convenient quantity 
for characterizing the perturbations is the intrinsic cur- 
vature perturbation of the comoving hypersurfaces |30| . 
£ = u/z, where u is a gauge-invariant scalar perturba- 
tion 0, and z ee a/(c s H)[-H/{AnG) +K,/a] 1 /'\ where 
c s is the sound speed and JC is the curvature of spatial 
sections. For single-field inflationary models, this simpli- 
fies to z = a(j)/H. The quantity u satisfies the dynamical 
equation 



A. Background Equations 

For single- field inflationary models, i.e., when inflation 
is driven by a single homogeneous scalar field, the dy- 
namical equation for the inflaton field </> is given by 



4>{t) + 3H(t)4>(t) 



dV(4>) 
d<f> 



= 0. 



(1) 



The dots denote derivatives with respect to physical 
time t. The evolution of the Hubble parameter H(t) = 
a(t)/a(t), where a(t) is the scale factor, is given by the 
Fricdmann equation 



8ttG 



~j?{t) + V{4>) 



(2) 



Equivalently, by taking the derivative of H(t) with re- 
spect to time and inserting the equation of motion Q 
for d> we can write 



Hit) = -4irG<j) 2 (t) 



(3) 



Here and in the following we have set H = c = 1. 
Eqns. QJ and 101 or © determine the evolution of the 
inflaton field completely. In addition, it is useful to cal- 
culate the number of expansion e-folds via 



N(t) = H(t). 



(4) 



Monitoring the number of e-folds for which inflation lasts 
allows us to select observationally relevant inflationary 
models and their parameters. Finally, the conformal time 
r](t) is defined by 



a(t) 



(5) 



Au - 



0, 



(6) 



where primes denote derivatives with respect to confor- 
mal time and A is the spatial Laplacian in comoving 
coordinates. It follows immediately that £ is approx- 
imately constant in the long wavelength limit k — > 0. 
This is true during the inflationary phase as well as dur- 
ing the post-reheating era. Moreover, the Einstein equa- 
tions can be used to connect the gravitational potential 
<&A and C so that a computation of the power spectrum 
of £ provides all the information needed (aside from the 
transfer functions) to extract the temperature anisotropy 
of the CMBR. Details of this procedure can be found in 
Refs. jMEJ. 

The calculation of the relevant power spectra involves 
a computation of the two-point functions for the appro- 
priate quantum operators, e.g., 



(0\u(r), x)u(r), x + r)|0) 

the operator u being written as 
d 3 k 



dk sin kr 
k kr 



Pu(v,k), (7) 



u{x) 



(2tt) 3 / 2 



a k u k (r,)e tk - x + aiulirit 



(8) 



where dk , d k are annihilation and creation operators with 
[a^,a|,,] — Skk', and dk\0) = Vfc. The complex ampli- 
tude Uk(r]) satisfies 



u k = 0. 



(9) 



Solving Eqn. is the fundamental problem in determin- 
ing the primordial power spectrum. The corresponding 
mode equation for tensor perturbations is given by 



As explained later, the uniform approximation is ex- 
pressed naturally in terms of the conformal time. 

Once the background equations are solved we can pro- 
ceed to the next step, the evaluation of the scalar and 
tensor perturbations. 



a 
a 



4 +{k 2 -—)v k = 0. 



(10) 



Once the mode equations are solved for different mo- 
menta k the power spectra for the scalar and tensor 



4 



modes are obtained via 



III. THE UNIFORM APPROXIMATION 



P s (k) = lim — 



P T {k)= lim ^ 

fcr;->0- 2-K 2 



z{ri) 



a{r)) 



(11) 
(12) 



where we denote the power spectrum for £ by Pg. 

The tensor power spectrum is often defined with an 
additional factor as 



P h {k) = 8P T (k), 



(13) 



leading to the definition of the tensor to scalar ratio as 
P h (k) 8P T (k) 



R{k) 



P s (k) P s (k) 



(14) 



The generalized spectral indices for the scalar and ten- 
sor perturbations are defined to be 



n s (k) 
n T {k) 



= 1 



dlnP s (k) 



dink 
d\nP T (k) 
dink 



(15) 
(16) 



Running of the spectral indices is conventionally parame- 
terized by the second logarithmic derivative of the power 
spectra: 



a s (k) 
ar(k) 



d\n ns(k) 
dink 

din nx{k) 
dink 



(17) 
(18) 



Numerical results for the power spectrum are often fit 
in the literature assuming a power-law behavior and a 
small running of the spectral index around a pivot scale 
k* (see e.g., Refs. [25| H2)- This fitting ansatz for the 
power spectrum is 



P(k) = A 



fit 



rifit + 5 "fit In - 



(19) 



the fitting parameters being Afn, «fit and aat- The spec- 
tral index is evaluated at the pivot scale: 



ns(k*) = 1 + nfi t . 



(20) 



The running is parameterized by agt. This fitting form is 
not self-consistent since it is not possible to have strictly 
constant n and non-zero a given the above definitions. 
This inconsistency manifests itself by the uncontrolled 
growth of errors away from the pivot scale. 

In the following we will describe how to obtain approx- 
imate solutions for the scalar and tensor power spectra 
and the spectral indices and their running, assuming that 
the evolution of the background quantities is known. 



In two recent papers 0, lldj| we have implemented an 
approach based on the uniform approximation (for a de- 
tailed description of the approximation, see Ref. [lj|) to 
calculate the power spectra and the associated spectral 
indices of primordial perturbations from inflation. Our 
method leads to simple expressions for the power spectra 
and spectral indices with calculable error bounds. The 
solution for the modes Uk can in principle be determined 
to arbitrary order in the uniform approximation. In prac- 
tice, results accurate to sub-percent level are obtained at 
next-to-leading order. 

While a detailed derivation of the mode equations, 
power spectra, and spectral indices with the correspond- 
ing error terms is given in Ref. |10|. for completeness 
and to set notation, we provide a brief summary of the 
main equations. In addition, we give an expression for 
the scalar to tensor power spectrum ratio R(k) with the 
corresponding error term and explain how to compute 
the error terms in the uniform approximation in an ef- 
ficient way. A new result, shown here, is that a simple 
improvement of the leading order uniform approximation 
leads to very good accuracy, generically better than 0.1%. 
This improvement is obtained by utilizing knowledge of 
the next-to-leading order results without actually imple- 
menting them fully. 



A. The Mode Functions 

In order to solve for the scalar and tensor mode func- 
tions in the uniform approximation, it is necessary to 
rewrite the differential equations 10 and l|10|) in the fol- 
lowing form: 



1 



4(v) 



1 



where v 2 s = (z"/z)f/ 2 + 1/4, and 



-k 2 + ^ 



'(»/) - 7 



(21) 



(22) 



where v\ 



{a"/a)rf + 1/4. The shift of 1/4 in the 
u>jiiiiiliuii ui iSg and is necessary in order to have a 
convergent error control function [lOj. In the following 
we will describe the analysis for the scalar modes and 
present only the final results for the tensor modes. 

The general solution for Uk(r]) is a linear combination 
of the two fundamental solutions u^\rj) and u^\rf), viz., 



u k {r l )=A{k)u^\r l )+B{k)u^{r l ), 



(2) 



(23) 



independent of the order of the approximation. To fix 
the coefficients A(k) and B(k), a linear combination 

of u^\r]) and (rj) must be taken so that Uk(i]) — 
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e tkri /y2k in the limit krj — > — oo. With proper normal- 
ization, the solution for is, in leading order, 

Ufc.i.sfa) - yf ^)ff5 V4 (^ »?) [Ai(/ $ ) - iBi(/ s )] , 

(24) 

with 

{o /.jjg i 2/3 

±2/ d V '[T9s(k,v')] 1/2 j ,(25) 

S5 (M) = #^-fc 2 . (26) 

One part of the solution is valid to the left of the turn- 
ing point fjs, defined as the solution to k 2 = Vg{fjs)/fjg, 
and the other part is valid to the right of the turning 
point (77 > fj). The uniform approximation allows us to 
calculate bounds on the errors. We write 



In a similar fashion the results for the tensor modes Vk 
can be derived. We have 



= fcJS- - lC \l 2^k) 6XP { 3 ^.>(*'^ >/a 

(32) 

where /t,> indicates that we replace gs{k 1 rj) with 
gr(k, rj) and fjs with 77T in Eqn. 1)25(1. 

Having obtained expressions for the scalar perturba- 
tions Uk and the tensor perturbations Vk we can now de- 
rive the corresponding scalar and tensor power spectra. 



B. The Power Spectra 

The expression for the scalar power spectrum Ps(k) as 
defined in Eqn. (TJ| in the uniform approximation with 
the corresponding error term is 



(27) 



where the error term encapsulates the contribution to 
Mfc.i. g b eyond leading order. As derived in detail in 
Ref. |10| , the error term is bounded by 



|e M , $ (r?)| < ^{[exp(AV^(f))-l] 
+ [exp(AV^(£))-l]}, 



(28) 



where V(£ ) denotes the total variation of the error control 
function £{rf). A numerical estimate shows A ~ 1.04 [Isj . 
The error control function reads 



J fjs Ussl 



1 



/4 <V 2 Vlssl 1 / 4 . 
1 5| 5S | 1/2 



4r?' 2 | 5s |V2 i6|/ SiS | 



d??'. (29) 



The krj — > 0~ limit defines the region of interest for 
calculating power spectra and the associated spectral in- 
dices. In this region, the I/77 2 pole dominates the behav- 
ior of the solutions and the Airy solution goes over to 
the LG form leading to simple expressions for the spec- 
tral indices 01 • The region of interest lies to the right of 
the turning point where the argument of the Airy func- 
tions becomes large. This allows the approximation of 
the Airy functions in terms of exponentials, leading to 



Uk,i,>{v) 



-^=5s 1/4 (M) 



2 6XP 



--[/ S ,>(M)] 3/2 



-jexp<{ I [fs,>(k,v)f /2 



(30) 



For computing the power spectra in the kr\ — ► limit, 
only the growing part of the solution is relevant: 



Uk.i,>(v) = l™ —iC 



-v 



2Mv) exp t|[ /s '>^ 3/;i 



(31) 



Ps(k) 



k 3 

lim — -r- 

fctj->o- 2tt z 



Uk,i,>(v) 



lim P hS (k) [1 



\l + e k ,i,>(v)\ 
k,i,sM] > (33) 



with 



JP 

e fc,l,S 



2Re e k ,1 ,> + |ejt,i>>l 



(34) 



where Pi,s(k) denotes the power spectrum for the scalar 
perturbations in the leading order approximation. We 
can now substitute either the full expression for Uk given 
in Eqn. or the LG form from Eqn. IpHJl. 
Using the LG expression for u k we have 



Pi,s{k)= lim 



1 



kr,^o- 4tt 2 I z (77) 1 2 1/5(17) 



exp 



[fs,>(k,r])] 



3/2 



(35) 

with the error term for the power spectrum given in 
Eqn. 134fl . The calculation for the tensor power spec- 
trum follows along the same lines, yielding 



P hT (k) = lim 



-V 



kn^o- 47r 2 |a(7y)| 2 vt(v) 



exp 



[/t,>(M)] 



3/2 



(36) 

with the error term in the same form as in Eqn. 134|) with 
the substitution S — > T. 

The tensor to scalar ratio, R(k), is given by 



R(k) = 

with the error term 



8Pi,r(fc) 
Pi,s(k) 



(1 



e k,l 



c k,l,T 



1 



-k,l) 



1. 



(37) 



(38) 



-fc,l,S 



In the case of power- law inflation, where v is constant, 
the error is identically zero, as in this case v$ — z/y. In 
other words, the ratio of tensor to scalar perturbations 
for power-law inflation is exact already at leading order 
in the uniform approximation. 
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C. Estimate of the Error Bound 

Although we can calculate the error bound for the 
power spectra from the general expressions in Eqns. i|28|) 
and (|29[1 . it is convenient to have simpler estimates for 
the errors. 

We begin by considering the case of constant v, where 
the fc-independent error bound for the power spectrum 
in leading order of the uniform approximation is [l^ 



A 



1 



> ' W ' 36^2 1 J ' (39) 

the generic v denoting either of vs or vt- This bound is 
rigorous and useful, though the prefactor is not optimally 
sharp for the case of constant v. 

Suppose now that v(r\) varies slowly with time. Fix 
k and consider the value of v(rf) at the turning point 
fj(k), defined by kfj = —i>(fj). Given the slow variation 
of ^ (r?), this value 9(k) is a slowly varying function of k. 
By expanding the expression for the error control func- 
tion of Eqn. H29|) locally around the turning point, we 
obtain what is in effect a derivative expansion for the 
error term. The leading term in this expansion, which 
is free of derivatives, has the same form as the expres- 
sion above for constant v, though it now carries the mild 
fc-dependence of the variable v case, 



c /c,l 



< 2V2 



1 



A 



+ 



69(k) 72P 2 (fc) 36V2P 2 (fc) 



This bound is not meant to be rigorous, since higher or- 
der terms in the derivative expansion are not included. 
However, it is effective and useful in the case of slowly 
varying v{r\). For further discussion of local approxima- 
tions of this type, see also Section lill E 21 In Ref. [lOj we 
have extensively discussed how a non-constant v gives 
rise to a non-vanishing and fc-dependent error for the 
spectral index. 



The natural question arises as to whether or not it is 
possible to utilize these results to improve the leading 
order expressions without recourse to full computation 
of the sub- leading approximations. For most viable in- 
flationary models, v varies slowly and corrections from 
the derivatives of v are sub-dominant. The full next-to- 
leading order machinery may not be required if v is suf- 
ficiently well-behaved. In this section we implement this 
idea and derive improved leading-order results for the 
power spectra. Note that, at leading order, the spectral 
index is exact for constant v, thus the main improvement 
to be expected is in the amplitude of the power spectrum. 

We begin with an argument from our previous 
work , which was used to understand why the leading- 
order result for the spectral index was much more precise 
than estimated bounds for the power spectrum such as 
given in Eqn. 14U|) would seem to indicate. The error is 
in fact dominated by an amplitude prefactor which has 
only a subdominant contribution to the spectral index. 
We split the scattering potential in the form 



- 2 (v) - 7 



(42) 



and choose the ^-independent but fc-dependent constant 
9(h) to be the value of 1/(77) at the turning point. As 
argued previously , this splitting allows us to identify 
two separate contributions to the total error term for the 
power spectrum, 



(43) 



The term e arises solely from the ultra-local contribution 
in the derivative expansion of v 1 and by explicit calcula- 
tion it is known to be of the form 



i = [r*(p)] 2 -i 
1 1 



(44) 



31 



139 



69(h) 72P 2 (fc) 6480i> 3 (fc) 155520P 4 (fc) 



D. Improvement of the Leading Order Result 

In this section we present a new improvement for 
the leading-order uniform approximation. In previous 
work [lfjj we derived next-to-leading order results for Uk 
and i>k and the corresponding power spectra and spec- 
tral indices. The resulting expressions contained several 
integrals, which are tedious to evaluate. On the other 
hand, the results at next-to-leading order for the case 
of constant v turned out to be very simple. In essence, 
higher order terms in the uniform approximation gen- 
erate a prefactor which occurs in the Stirling series for 
Y(v) 



a p 

m 



r» = i + — + — 

w 12i/ 288^ 



139 



(41) 



51840^ 3 

which we know a-priori to be present in the case of con- 
stant v. This all-orders prefactor improves the normal- 
ization of the power spectrum dramatically. 



The remaining error term satisfies an integral equation 
of the form considered by Olver [Tsj , with a reduced in- 
homogeneity; explicit calculation leads to the rigorous 
bound 



|g[< V %J ] [l + 0(\e\)]\e\ 



(45) 



where £ is the full error control function, £ is the error 
control function for the case of constant v = 9(k), and 
V(-) indicates total variation as before. 

In the generic case of slowly varying 1/(77) this expres- 
sion clearly shows that e is significantly reduced in com- 
parison to e. Therefore we are motivated to absorb the 
ultra-local contributions from higher order corrections 
into the power spectrum, leading to an improved first- 
order expression 



A(fc) = F 1 (fc)[r*(p)] 2 



(46) 
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This improvement applies to both the scalar and tensor 
power spectra, with vs{v) and vt{v) respectively. Since 
typical values of v ~ 2 occur for both scalar and tensor 
fluctuations, the terms on the right hand side of Eqn. i|46|) 
correspond to amplitude corrections of order 10%, 0.35%, 
0.06%, and 0.006%, respectively. These can be viewed 
primarily as local normalization corrections. 

To complete the discussion, we can obtain a non- 
rigorous estimate for the size of e, again using a deriva- 
tive expansion in £ to isolate the leading local contribu- 
tions. In this expansion, the leading derivative-free terms 
in £ — £ must cancel, and therefore the leading term is 
proportional to the first derivative of v(rj). After some 
calculation we find 



3 

< - 
~ 2 



< 



[1 + C(|e|)]|e 
d\n.v(k) 



AD{k) 



d\nk 



l + O 



(47) 



where we have used the chain rule to write the derivative 
in terms of a derivative of v(k) with respect to k and have 
used the explicit form for e(fc). It is understood that this 
is not a rigorous inequality, since we have neglected the 
higher order terms in the derivative expansion, but it is 
well-motivated and useful for slowly- varying v(ri). This 
error estimate determines the order to which the local 
corrections (I46|l should be taken into account. For in- 
stance, when v is slowly varying, most of the error can 
be compensated using Eqn. I|46l) . Examples will be en- 
countered in Section lVll where we explicitly demonstrate 
the success of this improvement procedure. As is to be 
expected, the ratio R(k) is much less sensitive to the nor- 
malization error compared to the amplitude of the power 
spectrum. For this quantity, the error is well-estimated 
by Eqn. (|?7|) . 



E. The Spectral Indices 

We now turn to the evaluation of spectral indices and 
corresponding error terms from the power spectra com- 
puted using the uniform approximation. We will discuss 
further simplifications leading to local expressions for ns 
and nx- 



1. Integral Expression for the Spectral Index 

The differentiation of the power spectrum with respect 
to k is straightforward. As stated earlier, it is important 
to remember that the turning point fjs is a function of 
k since k — — vs(jjs)/\vs\ where i>s(vs) is the value of 
1/5(77) at the turning point T) = fjs- Using this relation, 
we find 



n hS (k) = 4 -2k 2 



lim 

krj— >0- 



drf 



fjs V9s{k~rf) 



(48) 



The error for the spectral index only arises from the 
fc-dependent part of the error in the power spectrum. 
Therefore, the error in the spectral index is sensitive only 
to the time variation of vs- To estimate this error, the 
spectral index as written in Eqn. I|15f) can be expressed 
via the leading order power spectrum in the following 
form: 



n s {k) 



1 



dlnPig 



d\n k 



de ii,s 
' dk 



with 



t k,l,S 



de k.l.S 

3 dk : 



(49) 



(50) 



where we estimate 1 s by the leading order error term 
in Eqn. (|i2"|l . i.e., by e as given in Eqn. It should 

be noted that the leading errors in the spectral index are 
proportional to k derivatives of v{k), unlike the situation 
for the power spectrum [Eqn. 14411 ]. 

The above analysis can be carried out for tensor per- 
turbations in an identical fashion, including the error es- 
timation, with the replacement z/g — > vt- The spectral 
index for gravitational waves is given by 



ni )T (k) = 3 



2k 2 lim 

fcrj— >Q~ 



vt \/gr(k,r]') 



(51) 



2. The Local Approximation 

It is possible to simplify the expressions for the spectral 
indices as given in Eqns. (|48|l and (|51|l by approximating 
the integrals further. By doing so, we lose the ability to 
quantify the error estimates for ng and tit but are able 
to write down local expressions for the spectral indices. 
The integrands in Eqns. (|48() and (|51|l have square-root 
singularities at the turning points, i.e., at the lower inte- 
gration limits. At the upper limit, r\ goes to zero and the 
integrands vanish linearly, therefore, assuming v$ (r)) and 
vt(v)i respectively, are well-behaved, we expect the main 
contribution to the integrals to arise from the lower limit. 
Using the knowledge that vs and vt are slowly varying, 
we expand them around their turning points in Taylor 
series. To second order in derivatives, 

- v 2 s +2V s D' s (r? - fjs) + {v's +vps)iv-Vs) 2 i (52) 

where a bar indicates that the quantity has to be eval- 
uated at the turning point. For vt we obtain a similar 
expression with the index S replaced by T. We can now 
solve the integrals in Eqns. (|48|l and (|51|l exactly and find 
for the scalar spectral index 



n s (k) ~ 4-2P s \l 



-Vs 



(53) 



+ 



Vs 



'I 



I'Z 



(2-tt) + ^(1 



vs 
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and for the tensor spectral index analogously 
n T (k) ~ 3-2^ T |l-^ T (l-|) 



(54) 



2 



^(2-^ + ^(1-^) 



ISrp 



The local approximation is a further simplification of 
the leading order result in the uniform approximation 
and we will test its accuracy with three numerical exam- 
ples in Sect ion IVll We denote the leading contribution of 
unity in the curly brackets in Eqns. I|53ll and (|54(l as ze- 
roth order, the second term proportional to v' as first or- 
der in the local approximation, and the remaining terms 
proportional to v' 2 and v" as second order. 



Equivalently, for the tensor spectral index derived from a 
slow-roll expansion of the local result, Eqn. 154|) . we find 



/,,(/.:)~-2e- ( y . 



1 'f-Shr )<«>,.. 



where we have used 



a 
a 



2a 2 H 2 [ 1 - -e 



(60) 



Compared to results given previously 0, 0] , the expres- 
sions for the spectral indices given here have been fully 
expanded to second order in slow-roll parameters. 



IV. THE SLOW-ROLL APPROACH 



3. Spectral Indices in Terms of Slow-Roll Parameters 

It is possible to simplify the expressions for the spec- 
tral indices one step further and obtain results similar 
to the ones obtained in the slow-roll approximation to 
be discussed in the next section. We have inve stig ated 
this further approximation in detail previously | lf| and 
pointed out its shortcomings. Here the results are sum- 
marized for leading order only. 

We begin by considering the spectral index for the 
scalar perturbations as given in Eqn. JSHJ). In order to 
write ns in terms of slow-roll parameters which are de- 
fined as (see e.g., Ref. |34j 1 




1 d n+1 cj) 
1FZ dt^ 1 



(55) 



we expand v' s fjs, and Vgfjg up to second order in 
these parameters. The contribution v'gfjg/vs is already 
of fourth order and is neglected. Using the expression for 
conformal time in slow-roll parameters given by 



>1 



1 



3e 2 



2eS 1 



1 

"aH 

the relation, v\ = (z" / z)rj 2 + 1/4, and 

3 



2a 2 H 2 



2e£i 



(56) 



(57) 



it is easy to write down v and its derivatives in terms e, 8\ , 
and &2 (For a detailed derivation and explicit expressions, 
see Ref. ^3). As m earlier sections, the bar indicates that 
the parameters are to be calculated at the turning point. 
Inserting all the expressions into the local approximation 
for the scalar spectral index, Eqn. i|53|) , allows us to write 
the spectral index in terms of slow-roll parameters: 



ns(k) 



l-4e-2<5i 



(58) 



2(St - S 2 ) — - 7T 
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The most popular analytic method to evaluate primor- 
dial power spectra and spectral indices is the so-called 
slow-roll approach. The basic idea behind this approach 
is the following. Recall the expression for z" / z in terms 
of e, 5i and 82 as given in Eqn. 1)57(1. Assume now that 
e and 8± are small and constant. If this is the case, all 
terms 8 n> \ vanish, since they can be written as deriva- 
tives of e and 8\, and all terms of order 0(e 2 ) can be 
neglected. It is then also possible to express a 2 H 2 as 
(1 + 2e)/rj 2 , leading to a simplified equation for U}~: 



with 



4 + (k 2 - A/if) u k = 0, 



A = 2 + 6e + 3^i = constant. 



(61) 



(62) 



This equation is immediately solved in terms of Bessel 
functions. Since the constant A will be different for differ- 
ent values of k (as e and <5i are approximated by different 
constants for different k), for every momentum k a "best- 
fit power-law" result is obtained and expressions for the 
power spectrum Ps(k) and the spectral index ns(k) can 
be derived. Tensor perturbations can be treated in the 
same way, leading to expressions for Pr{k) and nxik). 
The terminology "slow-roll" is easy to understand: the 
requirement on e and 81 being small and roughly con- 
stant admits inflationary models with smooth potentials, 
leading to a phase where the inflaton field rolls slowly 
down the potential before inflation ends and reheating 
begins. The shortcoming of the slow-roll approach is 
that there is no simple way to improve the approxima- 
tion beyond leading order. Once the assumption that 
the slow-roll parameters are constant is given up, at the 
next formal order of the approximation, the Bessel solu- 
tion is no longer valid [35j . Different ways for improving 
the slow-rol l ap proximation have been suggested, e.g., 
[34l l36l I3H l38| . but they often lead to rather involved 
expressions for power spectra and spectral indices and, 
more importantly, are not error-controlled. 

As the slow- roll approach is convenient, widely used, 
and for most common inflationary models expected to 
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be a reasonable approximation, we will compare slow- 
roll results against the exact mode-by-mode integration 
and the uniform approximation in Section IVII For com- 
pleteness, we summarize here the results for the power 
spectra and spectral indices following closely the work of 
Stewart and Lyth [39j |. 

The power spectra in the slow-roll approximation read 



P s (k) ~ [l + 2(2-ln2-&)(2e + £i)-2e] 
2GH 4 



iH=k 



P T (k) ~ [1 - 2(ln2 + 6- l)e] 



2GH 2 



(63) 
(64) 



iH=k 



where b is the Euler-Mascheroni constant, 2 — In 2 — b 
0.7296 and In 2 + b - 1 ~ 0.2704. 

The scalar and tensor spectral indices are given by 



n s {k) ~ l-4e-25i-2(l + c)e 2 + i(3-5c)e(Si 



n T {k) 



-2e- (3 + c)e 2 - (1 + c)eS u 



(65) 
(66) 



and c ~ 0.08145. The expressions ns(k) ~ 1 — 4e — 28 1 
and rixik) ~ — 2e are the well-known first-order slow-roll 
results. 

In the numerical calculations for slow-roll presented 
below, we will approximate the solutions of the mode 
equations in terms of Bessel functions, while solving the 
background equations without any further approxima- 
tions. Hence, we do not assume the strict slow-roll con- 
ditions that the slow roll parameters e, 8±, and 82 be 
constant. This is not the way the slow-roll approach 
is often implemented when estimating inflationary pa- 
rameters: A further asymptotic approximation is utilized 
to solve the background equations 0> EH employing a 
derivative expansion in the inflationary potential V (</)). 
Such expansions may provide useful estimates but do 
not give the precision required by present and upcom- 
ing CMBR observations. The reason for the inaccuracy 
arising from expansions in the potential, as will be seen 
in an example, is that the Taylor expansion in derivatives 
can often fail. This is not the case when the parameters 
are expressed in terms of the Hubble parameter and its 
derivatives. Our aim here is to determine the accuracy of 
the mode solutions themselves, therefore we do not focus 
on errors arising from approximations to the background 
equations. 



V. NUMERICAL IMPLEMENTATIONS 

In this section we describe the numerical implemen- 
tation of the exact mode-by-mode integration and the 
uniform approximation and its simplifications. Technical 
details can be found in Appendix^ Since the uniform 
approximation is useful both in setting initial conditions 



for the exact numerical solution, and as a complete semi- 
analytic method in its own right, we discuss it first below. 



Leading Order Uniform Approximation: 
Numerical Issues 



1. Preliminaries 

We begin by addressing some technical points and in- 
troducing conventions. Since the conformal time is de- 
fined only up to a constant, we set it to zero at the end 
of inflation. If the given model does not have a natural 
end to inflation this point is somewhat arbitrary and we 
set i] to zero typically some number of e-folds after the 
highest mode of interest freezes out. 

The power spectrum and the spectral indices are to be 
calculated in the limit kr\ — > 0~, so that the LG approxi- 
mation can be used. To avoid accumulation of numerical 
error, however, these quantities should not be calculated 
directly in this limit. Since the freeze-out happens soon 
after the turning point is crossed, the computation is 
carried out some e-folds after the turning point for each 
respective mode, but well before the end of inflation. The 
linear combination of both solutions to the mode equa- 
tions, in terms of Ai- and Bi-functions, is then completely 
dominated by the exponentially growing part. We found 
that carrying out the computations 4 or 5 e-folds after 
the respective turning point provided sufficient accuracy. 

In all the computations described below, the Ai- and 
Bi-functions were calculated with the algorithm given in 
Ref. HJ. 



2. The Mode Solutions 

The full solutions Uk,i,>(j]) an d u k,i,<(v) are n °t 
needed for the calculation of the power spectrum and 
the spectral index. Nevertheless, it is useful to calculate 
some of them for selected momenta k in order to compare 
the leading order uniform solutions to the exact numeri- 
cal solutions found from a mode-by-mode integration of 
the differential equations. For 77 < fj we can calculate the 
integrals appearing in f<(k,rj), as defined in Eqn. it2*5|) 
numerically via 



dr)'Vg{k^f)=( + Jdi/VKMO, (67) 



where r\i is an initial value of the conformal time. In 
the actual numerical routine we therefore have to know 
the second integral before we can calculate the uniform 
solutions left of the turning point: This is achieved by 
an additional run of the integrator for the background 
equations. 
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3. Power Spectra and Spectral Indices 

The power spectra for scalar and tensor perturbations 
in leading order of the uniform approximation are given 
by Eqns. 1(35(1 and ((36(1 . The integrals are calculated us- 
ing a trapezoidal rule with non-equidistant discretization 
in conformal time. As mentioned above, we avoid calcu- 
lating the spectra numerically in the limit krj — > _ , but 
instead do so some 4-5 e- folds after the turning point. 
In the case of power-law inflation analytical results are 
available for the leading order contributions to the power 
spectra (see Ref. ^3)- We have checked that the power 
spectra numerically calculated from Eqns. 1)35(1 and 1(36(1 
are in agreement with these analytic results (see Sec- 
tion EO] and Table H| for details). 

The spectral indices in the uniform approximation may 
be calculated either by numerically differentiating the 
power spectrum as described in Section IlII El for the ex- 
act mode-by-mode approach, or by using the formulae in 
Eqns. and (|S"T)l . In Appendix IA 2l we describe how to 
deal with the inverse square root singularities appearing 
in the integrals to be performed in the second case. 



B. Mode-by-Mode Numerical Integration 



0.00 
|« -0.05 

s 

2 -o.io 

I -°- 15 

^ -0.20 





Analytical result 




Numerical result 


■ 


Uniform, 1. Order 







-4 -2 2 4 

e-folds 

FIG. 2: Imaginary (growing) part of the scalar mode func- 
tion for a power- law model with p = 11 for the mode 
k — 0.0214 h Mpc -1 . The solid red line is the exact ana- 
lytical solution, the dashed-dotted blue line the result from 
the numerical solution of the exact equations, and the dashed 
green line the uniform approximation to leading order; the 
analytical and uniform approximation results are almost on 
top of each other. The arrow shows the turning point and the 
green band is the estimated error bound for the leading order 
uniform approximation. 



1. Initial Conditions and Mode Functions 

In order to numerically obtain the mode functions, we 
must satisfy the initial condition requirement, i.e., in the 
limit kr/ —i- — oo, 



Uk(v) 



1 



-ikr) 



2k 



(68) 



Two difficulties in imposing this formal initial condi- 
tion immediately arise. First, in any numerical solution, 
the calculation must begin at a finite initial time, thus 
for modes with small enough values of fc, the condition 
krj — > — oo may not be fulfilled. Second, for modes at 
larger k values, there are very many oscillations before 
the turning point is reached, naively requiring very fine 
time-steps if the entire temporal range must be handled 
numerically. 

To circumvent these problems, we use the uniform ap- 
proximation to set initial conditions in a regime where 
it is exponentially accurate. For each mode, we take 
as initial condition the uniform approximation result at 
roughly 20 zeros, i.e., 10 oscillations, before the turn- 
ing point for that mode. The number of zeros from 
a given time to the turning point can be estimated by 
nix w k[f)(k) — rj\. As shown below, this initialization pro- 
cedure suppresses numerical errors, especially in the high 
k regime where the precision of the power spectrum and 
the spectral index is improved without taking smaller and 
smaller time steps. In addition, only a smaller number of 
"active" modes, i.e., the modes that are within some 20 
zeros before the turning point and not yet frozen, need be 
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FIG. 3: The absolute value of the scalar mode function for 
power-law model following Fig. |21 



considered at any time. In the regime of small fc, we first 
note that as inflation has to start somewhere in practical 
numerical calculations, this introduces a lowest value for 
k, defined by the criterion that the mode should be well 
inside the Hubble length. As the initial conditions for 
inflation are unknown, we will assume (i) that inflation 
began well before the 55-65 e-folds necessary to solve the 
flatness and horizon problems, and (ii) by the time our 
calculations are to be performed, Eqn. ((68() applies. By 
isolating how initial conditions are defined from possible 
early-time artifacts, our method of implementing initial 
conditions also leads to substantial improvement in the 
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low-fc regime as well. 

In Figs. El and 01 we display the imaginary (grow- 
ing) part and absolute value of Uk for a mode with 
k = 0.0214 h Mpc -1 for a power-law inflation model. As 
can be seen in both cases, the uniform approximation at 
leading order and the numerical results are clearly very 
close. The error in the absolute value \uk\ (relevant to 
the power spectrum) from the mode-by-mode integra- 
tion is shown in Fig. The numerical error is less than 
1 part in 10 5 . To check the accuracy of the solution of 
the background equations, the numerical deviations from 
the expected pure constant values for v, e, and 5\ are 
shown in Fig. the errors are comfortably below parts 
per million. Detailed quantitative results are given in the 
following section. 




FIG. 4: Relative error of \uh\ shown in Fig.^lfor the numerical 
calculation (blue line). The arrow shows the turning point for 
the mode. 



0.0001 



le-05 



le-06 



le-07 




FIG. 5: Relative errors of v, e, and Si; model parameters are 
as specified in Fig. [5] 



2. Power Spectra and Spectral Indices 

For the case of scalar perturbations the modes freeze 
out once the power spectrum Pg(k,r]), as defined in 
Eqn. <|ll|l . becomes constant as a function of conformal 
time, i.e., when P' s (k,r]{ reezG ) ~ at a numerically de- 
termined freeze-out time 77f re eze- We found that it is nu- 
merically robust to track the time derivative of Ps(k,rj) 
to determine this freeze-out. The same situation holds 
for the tensor perturbations and Pr(k,r]) as defined in 
Eqn. l(T2l . 

Once the power spectra have been obtained, the spec- 
tral indices for scalar and tensor perturbations are found 
by evaluating the derivatives of the power spectra with 
respect to fc, as defined in Eqns. (|15|) and Ijlbj) . The 
derivatives are computed numerically with nonequidis- 
tant momentum discretization due to the momentum 
readjustment described in Appendix I A II We take three 
discretization points and approximate the derivative by 
a non-symmetric, second-order, rule. In the numerical 
evaluation of the derivatives, for every individual k we 
use two very close neighboring points with A In k sa 0.01- 
0.02 over the entire fc-range considered. 



VI. RESULTS 

In this section we discuss three inflation models. We 
start with power-law inflation, one of the very few ana- 
lytically solvable cases. Next we investigate two chaotic 
inflation models, with a quadratic and quartic potential, 
respectively. Finally, we present a model in which the 
third derivative in the potential V(cf>) does not exist at 
one point, leading to a "kink" in z" / z. 

We present below results from numerical solution of 
the mode equations, the uniform approximation with dif- 
ferent simplifications and improvements, and slow-roll 
approximations in first and second order. In all cases, 
we show the results for the ratio of tensor and scalar 
power spectra, as well as the scalar and tensor spectral 
indices. Some of the results are presented in tables, for a 
specific value of the momentum fc, while selected results 
are shown in figures. In addition to the basic quantities 
themselves, we also give relative errors in some cases. 



A. Power-Law Inflation 

As one of the few analytically tractable models, power- 
law inflation |44| provides a useful foil for testing approxi- 
mations. This feature has maintained its popularity, even 
though the basic model is not realistic, as inflation never 
comes to an end. We compared detailed results from the 
uniform approximation in leading and next-to-leading or- 
der with slow- roll and exact results previously The 
spectral index was exact already in leading order in the 
uniform approximation, while the slow-roll approxima- 
tion yields a Taylor expansion of the spectral index. The 
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TABLE I: Numerical precision tests in the power-law case with p = 11 at fc* = 0.06875ft Mpc" 1 = 0.0495 Mpc" 1 (h = 0.72±0.05; 
the WMAP pivot scale is at 0.05 Mpc -1 , see Ref. Il5f i in the various approximations. 



Approximation 


R(k t ) 


ns(k*) 


«s(fc*) 


n T (fc,) 


a T (k*) 


Analytical 


1.454545 


0.8 





-0.2 





Numerical" 


1.454544 


0.79998 


< 10 7 


-0.20002 


< 10 7 


Uniform, 1. order 6 


1.454541 


0.79992 


0.00001 


-0.20008 


0.00001 


Local approx., 0. order 6 




0.80000 


< 10" 7 


-0.20000 


< 10" 7 


Local approx., 1. order 6 




0.79999 


< 10~ 6 


-0.20001 


< HT 6 


Uniform, Slow-Roll redux c 




0.80165 


< 10" 7 


-0.19835 


< 10" 7 


Slow-Roll, 1. order d 


1.454545 


0.81818 


< 10~ 7 


-0.18182 


< 10~ 7 


Slow-Roll, 2. order 6 




0.80165 


< 10" 7 


-0.19835 


< 10" 7 



a mode-by-mode integration 
'already exact in that order 
c see Eqn. and Eqn. S5§1 

d expected results: ng = 1 — 2/p and = —2/p 

e expected results: n$ = 1 — 2/p — 2/p 2 and = —2/p — 2/p 2 



amplitude of the power spectrum in the uniform approx- 
imation was accurate to roughly 10% in leading order 
while in next-to-leading order the accuracy of the ampli- 
tude improved dramatically with an error smaller than 
1%. These results hold for widely differing values of p, 
the power in the time evolution of the expansion a cx t p 
[the potential is V(cj>) — Vq exp(-\/2/p (/>)]. In contrast, 
the slow- roll approximation being an expansion in l/p, 
has large errors for small p, decreasing as p increases. 

Here, we use the power-law model for testing the accu- 
racy of our numerical implementations of the exact mode 
equations, as well as the uniform and slow-roll approxi- 
mations. The spectral indices are constant, therefore the 
running of the spectral index - which measures the k- 
dependence of ns and tit - vanishes. In addition, the 
ratio R(k) — 16/p of the power spectra is constant. In- 
stead of showing our results graphically - all curves would 
be almost indistinguishable by eye - we summarize our 
findings in Table U] We quote the results for R(k), the 
spectral indices, and their running as and ax at the 
WMAP pivot scale fc» = 0.0495 Mpc" 1 45]. We have 
picked the power p = 11, which is of course much too 
small to be considered a realistic inflationary model, but 
is perfectly acceptable for illustrative purposes. In the 
uniform approximation there is an explicit expression for 
the spectral indices ns and tt-t, while the tensor to scalar 
ratio R(k) is calculated via the amplitudes of the power 
spectra. In the purely numerical computation, the spec- 
tral index is obtained by taking numerical derivatives of 
the power spectrum as as described earlier. 

The error in the exact numerical determination of the 
ratio of the power spectra R*{k) is smaller than 1 part in 
10 6 , while the numerical error in the implementation of 
the uniform approximation is slightly larger (3 parts in 
10 6 ). In the case of power-law inflation, the spectral in- 
dices do not accurately reproduce the expected analytic 
results, though they do reproduce the expected slow-roll 



results to 1/p-order with the slow-roll parameters be- 
ing constant in time (in this case, e = —8% = l/p and 
62 = 2/p 2 ). This provides a check on the integration 
of the background equations, but has no further useful 
information. 

For the scalar spectral index ns(fc„) we find a rela- 
tive error in the exact numerical calculation of 0.0025% 
while for the uniform approximation the relative error 
is roughly 0.01%. Again, no error in the numerical im- 
plementation of the slow-roll approximation is detected. 
The situation for the tensor spectral index is sim- 

ilar. The numerical errors in the uniform approximation 
for the spectral indices are slightly larger than for the 
exact numerical mode-by-mode approach due to the nu- 
merical integration of a function with an inverse square 
root singularity. The results for as(fc*) and ar(k^) are 
very close to zero in all cases. 

In addition, we tested our improvement strategy for 
the uniform approximation which should work perfectly 
in the case of power-law inflation since v is constant. We 
found a relative error for Ps (k) of 0. 15% if we include cor- 
rections up to 1/z/ 2 , of 0.025% including corrections up to 
1/j/ 3 , and of 0.01% if we include corrections up to \/v A 
following Eqn. (|46[) . The corrections are slightly worse 
than anticipated but consistent with the small 0.01% er- 
ror in the numerical implementation of the uniform ap- 
proximation. 



B. Chaotic Inflation 

Next we study two chaotic inflation models. In these 
examples, the slow-roll approach is expected to work well; 
in order to obtain enough e-folds of inflation the poten- 
tials cannot be very steep, thus slow-roll turns out to 
be a very good approximation. For the uniform approx- 
imation we expect similar results as for the power-law 
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TABLE II: Determination of the characteristic quantities for the ^m 2 (f> 2 chaotic inflation model at fc* = 0.06875 ZiMpc 1 = 
0.0495±0.0034 Mpc -1 (h = 0.72±0.05; the WMAP pivot scale is at k = 0.05 Mpc -1 , see Ref 0) in the various approximations; 
parameters: m 2 = (1.89 ± 0.21) x 10" 12 /8ttG, (f)(0) = m.S/V^G, <f>(0) = -Q.I/VWkG/s leading to 57.655 e-folds of inflation 
after horizon crossing of fc=0.0495 Mpc -1 . 



Approximation 


R(k„) 


n s (k«) 


a s (k«) 


n T (k t ) 


a T (k t ) 


Numerical" 


0.13749 


0.96507 


-0.00064 


-0.01765 


-0.01779 


Uniform, 1. order 


0.13740(14) 


0.96505(3) 


-0.00064 


-0.01768(1) 


-0.01773 


Uniform, improved 1. order 6 


0.13749(0) 


0.96507(0) 


-0.00064 


-0.01765(0) 


-0.01773 


Local approx., 0. order c 




0.96465 


-0.00066 


-0.01786 


-0.01799 


Local approx., 1. order c 




0.96501 


-0.00064 


-0.01761 


-0.01770 


Uniform, Slow-Roll redux d 




0.96566 


-0.00062 


-0.01757 


-0.01768 


Slow-Roll, 1. order e 


0.13752 


0.96523 


-0.00063 


-0.01741 


-0.01755 


Slow-Roll, 2. order" 




0.96507 


-0.00064 


-0.01764 


-0.01779 



a mode-by-mode integration 

'the improvement here is to second order in powers of 
Cf. Eqn. 
c see Eqn. 1531 and 1541 
d see Eqn. 1581 and Eqn . 1591 
e see Eqn. lool and 1661 




FIG. 6: (a) Scalar power spectrum Ps(k) and (b) tensor power spectrum Pr(k) for a chaotic |m 2 2 -model, parameters as 
specified in Table [H] Solid red line: exact numerical results, dashed green line: uniform approximation, dashed-dotted blue 
line: slow-roll; dotted purple line: second-order improved uniform approximation [Cf. Eqn. I0OJ]; the green band is the estimate 
for the error bound for the uniform approximation, Eqn. (14011 . The exact results and the results from the improved uniform 
and slow-roll approximation are visually on top of each other. The error estimate for the improved uniform approximation, 
Eqn. is so small as to be indistinguishable from the result itself. 



case: an accuracy more or less independent of the chosen 
parameters. 



1. Quadratic Potential: V((f>) = m 2 (f) 2 /2 

The simplest chaotic model is that for a free scalar field 
with mass m. In this model, normalizing the amplitude of 
the scalar power spectrum to the WMAP fit, as described 
in Appendix IA 31 is equivalent to fixing m 2 = (1.89 ± 
0.21) x 1Q- 12 /8ttG. The initial value for the inflaton 



field is given by (f)(0) = 16.8/V87rG with a small initial 
velocity of <p(0) = — 0.1/vo7rG. These parameters ensure 
that the inflationary phase lasts long enough to provide 
a realistic model, i.e., it leads to 57.655 e-folds after the 
k = 0.0495Mpc _1 mode crosses the Hubble length. 

First, we analyze the results for the power spectra 
Ps(k) and Pr{k), plotted in Figs. & and 03. The ex- 
act numerical results for the power spectra as defined in 
Eqns. I|ll|l and (|12f) are shown in red (solid line), the re- 
sults from the leading order uniform approximation given 
in Eqns. 135|) and (|36|l are shown in green (dashed line), 
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FIG. 7: (a) Scalar spectral index ns(fc) and (b) tensor spectral index nr(k) for a chaotic ^m 2 (f) 2 -model, parameters as specified 
in Table ITT1 Solid red line: exact numerical results, dashed green line: uniform approximation, dashed-dotted blue line: slow- 
roll; the green band is the error estimate for the uniform approximation. Unlike the case for the power spectrum, accurate 
results for the uniform approximation are obtained without recourse to the improvement procedure specified by Eqn. 1461 . 




FIG. 8: Tensor to scalar ratio R{k) for a chaotic ^m?(j> 2 - 
model; parameters as specified in Table [n] All three lines lie 
practically on top of each other; see Fig. [0] for the relative 
errors. 



the first order slow-roll results as defined in Eqns. I|63|) 
and l|t)4|) are shown in blue (dashed-dotted line), and the 
improved uniform approximation results [to second-order 
in the sense of Eqn. I|4tj|) for the scalar perturbations] 
are shown in purple (dotted line). The light green band 
displays our estimate of the error bound for the leading 
order uniform approximation as given in Eqn. <|40[) . The 
error estimate for the improved uniform approximation 
in leading order [see Eqn. l(4"?|) ] is so close to the result 
itself that the error band is not visible in this plot. We 
will use these color and linestyle assignments for the re- 
mainder of the paper. 



FIG. 9: Relative error of the tensor to scalar ratio R(k) for 
a chaotic ^m 2 cj) 2 -model; parameters as specified in Table ITT1 
Here the dotted purple line denotes the (second-order) im- 
proved leading order of the uniform approximation; the green 
error band is the error estimate for the first order uniform 
approximation, the darker green band is the error estimate 
for the improved leading order result. 



The exact numerical results, the improved leading or- 
der uniform approximation, and the slow-roll approxima- 
tion are almost indistinguishable by eye in Fig. The 
leading order uniform approximation deviates in the am- 
plitude from the exact numerical result by ~ 10% as 
expected. The improvement strategy performs just as 
predicted by Eqn. (|46|l : the second-order correction re- 
duces the error to ~ 0.1% while the fourth-order reduces 
it further to ~ 0.01%, for both scalar and tensor pertur- 
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TABLE III: Determination of the characteristic quantities for the j\(f> 4 chaotic inflation model at k* — 0.06875 ZiMpc 1 = 
0.0495 ± 0.0034 Mpc" 1 (h = 0.72 ± 0.05; the WMAP pivot scale is at 0.05 Mpc -1 , see Ref 0) in the various approximations; 
parameters: A = (1.75 ± 0.19) x 10 -13 , (f)(0) = 24/v%rG, 0(0) = -1/V8ttG/s, leading to 60.579 e-folds of inflation horizon 
crossing of fc=0.05/Mpc. 



Approximation 


R(k«) 


ns(fe*) 


cts(k„) 


n T (fc„) 


ar{k t ) 


Numerical" 


0.25963 


0.94999 


-0.00090 


-0.03356 


-0.01705 


Uniform, 1. order 


0.25948(113) 


0.94990(3) 


-0.00089 


-0.03367(2) 


-0.01694 


Uniform, improved 1. order 6 


0.25964(0) 


0.94999(0) 


-0.00089 


-0.03356(0) 


-0.01703 


Local approx., 0. order c 




0.94942 


-0.00092 


-0.03395 


-0.01726 


Local approx., 1. order c 




0.94991 


-0.00090 


-0.03340 


-0.01697 


Uniform, Slow-Roll redux d 




0.95081 


-0.00087 


-0.03368 


-0.01706 


Slow-Roll, 1. order 13 


0.25969 


0.95077 


-0.00087 


-0.03285 


-0.01669 


Slow-Roll, 2. order" 




0.95001 


-0.00089 


-0.03354 


-0.01703 



a mode-by-mode integration 

'the improvement here is to second order in powers of 1/u, 
Cf. Eqn. i46) 
c see EamTSSl and |PJ 
d see Eqn. 1581 and Eqn . 1591 
e see Eqn. loot and 1661 



bations. This consistently good behavior is due to the 
fact that v is varying slowly. Note that the improve- 
ment strategy should be carried out to roughly match 
the error estimate given by Eqn. I|47|) - fourth-order in 
this instance; beyond this point the error is dominated 
by other contributions. 

Next we investigate the scalar and tensor spectral in- 
dices as functions of k. The results are displayed in 
Figs.[7^ and[7|3. The shaded band represents the error 
estimate for the leading order of the uniform approxi- 
mation [see Eqn. (|50|l . calculated with the estimate in 
Eqn. 144|) ]. The error estimate for the scalar spectral 
index is of the order ~ 0.002%, therefore of the same 
order as the discrepancy to the numerical result. Note 
that the improved scalar and tensor spectral index calcu- 
lations, following the discussion in Section TlIIE II agree 
completely with the numerical results in Tabic [n] The 
deviation of the slow-roll approximation for the scalar 
spectral index in first order from the exact numerical re- 
sult is also very small, roughly ~ 0.02%. As stated before 
this is not surprising: the slow-roll approximation is ex- 
pected to work well for this type of model, where the 
slow-roll parameters e and Si are very small and almost 
constant. However, the tensor spectral index has an error 
of more than 1%. 

In Fig. [S] the ratio R(k) of tensor to scalar perturba- 
tions as defined in Eqn. 114|) is shown. The correspond- 
ing relative errors of the uniform approximation, of the 
slow-roll approximation, and of the improved uniform ap- 
proximation arc presented in Fig. |5J In Fig. [S] the lines 
representing the different approximations are practically 
indistinguishable by eye. In contrast to the relative er- 
rors in the power spectra (~ 10% without improvement 
and ~ 0.01% with fourth-order improvement), the rel- 
ative error in the tensor to scalar ratio R(k) is already 



smaller than ~ 0.07% (dashed line in Fig. ^ without the 
improvement, and smaller than ~ 0.001%, i.e., practi- 
cally exact, with second-order improvement (solid line in 
Fig.© (Cf. Section irrrnj l. There is no point m improving 
further since the result has already reached the estimated 
error threshold i|47|) beyond which improvement becomes 
incomplete. This behavior is a general feature of the ap- 
proximation when the inflation model leads to a slowly 
varying v(rj). 

Finally, we have listed the ratio of tensor to scalar per- 
turbations, the scalar and spectral spectral index, and 
their running in the various approximations or simplified 
approximations (introduced in Section llll E|) , in TablelTTl 
All quantities are again evaluated at k* — 0.0495 Mpc -1 . 
Error estimates for the uniform approximation are indi- 
cated by the numbers in brackets. For chaotic inflation 
the running of the spectral indices, ats(k) and ax(k) is 
non-zero. Note that the running of the tensor spectral 
index is roughly an order of magnitude bigger than the 
running of the scalar spectral index. 



2. Quartic Potential: V(cj>) = X(j} 4 /4: 

In the quartic model - relative to the quadratic po- 
tential - higher derivatives of the potential exist. Con- 
sequently, the slow-roll results in this case are expected 
to have a bigger error. In comparison to the quadratic 
case, for the model considered below, the errors are worse 
for the scalar spectral index (~ 0.02% versus ~ 0.1%) 
and comparable for the case of the tensor spectral in- 
dex (~ 2%). In contrast, the uniform approximation still 
provides an accuracy of a fraction of a percent. Although 
the slow-roll expansion for this model can be improved 
to second-order with notably better results, this behavior 
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FIG. 10: (a) Scalar spectral index ns(k) and (b) tensor spectral index nr{k) for the quartic potential potential X4> 4 , parameters 
as specified in Table ITTT1 Solid red line: exact numerical results, dashed green line: uniform approximation, dashed-dotted blue 
line: slow-roll; the green band is the error estimate for the uniform approximation. 



exhibits the general tendency of the slow-roll expansion 
if terms arising from higher order derivatives of the po- 
tential are significant. In terms of observational viabil- 
ity, the </> 4 -model is under strong pressure from combined 
analysis of the WMAP CMBR anisotropy data and data 
from galaxy clustering (see e.g., Refs. |4a IItI Esj ) . 

The results for the spectral indices are displayed in 
Figs. IlUb and HUb . As in the previous example the 
leading-order uniform approximation is very close to the 
exact numerical results. However, the first-order slow-roll 
result does not match as closely as for the quadratic po- 
tential. As done earlier for the </> 2 -model, we have listed 
the various characteristic quantities in Table UTTl 

C. Inflationary Model with a C 2 -Potential Function 

As the last example we investigate a toy model with a 
continuous potential function with continuous first and 
second derivatives and a jump in the third derivative 
at a specific value of <p. Generally speaking, dynamical 
changes in the potential of the field driving inflation can 
be induced by couplings to other degrees of freedom. For 
example, in hybrid models a phase transition is used to 
terminate inflation. If a dynamical transition happens at 
cosmologically relevant scales, i.e., well before the end of 
inflation, it leaves a clear signature in the power spectra 
and spectral indices. Such a transition may be natu- 
rally realized in multi-field models of inflation (see e.g., 
Ref. Other examples are steps in the potential (see 

e.g., Ref. [IE IH|)> leading to oscillations in the primor- 
dial power spectra and spectral indices. In such models 
v 2 cannot be considered as a constant, but can display 
sudden changes. 

Rather than taking one of the potentials mentioned 
above we consider here a toy potential that is smoother 



in the sense that oscillations in v 2 or z" / z are avoided: 

v>0) = i™ 2 ^ (a _i) + ^ TO 2^ (1 _ a)0 

+ iam 2 2 + _Lm 2 (l-a)0 4 , ( 69 ) 

V<(4>) = \m 2 <j> 2 , (70) 
where V((j>) = V>(</>) for > 0* and V{4>) = V<(cj)) 







v s 2 (u) 








z'Vtza 2 ) . 


1 1 1 1 



5 10 15 20 25 

e-folds 

FIG. 11: z"/(za 2 ) and v% and a" /a 3 and u T for the C 2 - 
potential; the point 4>* is reached at N ~ 12.4. The beginning 
in time of the numerical calculation is at N — 0. The infla- 
tionary attractor is reached at N ~ 1.8, checked by varying 
4>(0) and determining at which e-fold the cj> behavior becomes 
independent of the initial velocities. The quantities a" /a 3 and 
vt, relevant for tensor perturbations, are much smoother than 
the corresponding quantities z" / za 2 and us for the scalar per- 
turbations, leading to smaller errors in the approximations as 
discussed in the text. 
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FIG. 12: (a) Scalar power spectrum Ps(k) and (b) tensor power spectrum Pr{k) for the C 2 -potential in Eqns. 1691 and 1701 : 
parameters: a = -100, m 2 = (1.90 ± 0.21) x 10~ 12 /8ttG, 0, = 15.2/V&ivG, 0(0) = 17.5/v^rS, 0(0) = -0.2/V8^Gs. Solid 
red line: exact numerical results, dashed green line: uniform approximation, dashed-dotted blue line: slow-roll; the green band 
is the estimate for the error bound for the (unimproved) uniform approximation. Again, the exact results and the results from 
the improved uniform and slow-roll approximation are on top of each other. 
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FIG. 13: (a) Relative errors for the scalar power spectrum Ps(k) and (b) tensor power spectrum Pr{k) for the C 2 -potential, 
Eqns. lb9H and I7UI . In both cases, the light green band denotes the estimated error bound for the first-order uniform 
approximation I40H and the dark green band denotes the best estimated error for the improvement procedure 146H . The results 
are nicely consistent with these estimates showing where the second-order improvement can be enhanced by going to higher 
order, and where it cannot. 



for < 0*. The potential is constructed in such a way 
that V>(0*) = V<(0*), V£(0«) = 1^(0* ) and V^'(0*) = 
V£(4>*) but V"'{4>*) + V£'(4>*). Thus there is a finite 
jump in the third derivative of the potential. 



We present numerical results with parameters chosen 
specifically to demonstrate the general effect of a more 
rapidly changing v. During the evolution in this poten- 
tial the parameters e and Si are not constant (not even 
approximately); 8i cannot be considered small when the 
inflaton field fulfills > 0* (with the parameters below, 
\6i\ can be as large as 0.14). 



The parameters chosen to specify the model are: a = 
-100, m 2 = (1.90±0.21) x 10- 12 /8ttG, 0* = 15.2 /V&nG, 
0(0) = 17.5/v / 8ttG, 0(0) = -0.2/v^ttG/s. With this 
choice of parameters the number of e-folds counted from 
K = 0.0495 Mpc' 1 is 57.320. 

The relevant time-dependent terms and z" /(za 2 ) in 
the scalar and Vj. and a" /a 3 in the tensor mode equations 
are displayed in Fig. llll as a function of expansion e-folds 
(z" I z and a" /a have been divided out by a 2 to filter out 
the exponential growth of the scale factor). The point 0* 
is reached at N w 12.4 (note that in this plot N = de- 
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FIG. 14: (a) Scalar spectral index ns(fe) and (b) tensor spectral index riT(k) for the C 2 -potential in Eqns. 1691 and I70H in the 
region around the kink; parameters specified as in Fig. 1121 Solid red line: exact numerical results, dashed green line: uniform 
approximation, dashed-dotted blue line: slow-roll; the green band is the error estimate for the uniform approximation. 



fines the beginning of the numerical calculation). Due to 
the jump in the third derivative both quantities for the 
scalar perturbations display a kink at this point. The 
qualitative behavior is also different on either side of the 
kink. While z" / z has a kink, a" /a is completely well- 
behaved. Note that z"/z, as e.g., expressed in Eqn. i(57|) 
as an exact expression in terms of the slow-roll param- 
eters, is more sensitive to higher derivatives of the po- 
tential than a" /a [Cf. Eqn. (|6()(l ]. Thus we can expect 
the effects of the change in the potential at <f> = 0* to 
be amplified in the scalar power spectrum relative to the 
tensor power spectrum. 

The results for the scalar and tensor power spectrum 
are displayed in Figs. 112b and 112b. With the same con- 
ventions as in Fig. [ft] the different approximations (lead- 
ing and improved leading order of the uniform approx- 
imation and the slow-roll approximation) are compared 
to the exact numerical results. The scalar power spec- 
trum (see Fig. 112b ) shows a significant deviation from a 
power-law shape. Up to k ~ 0.025 /iMpc -1 the spectrum 
rises, reaches a maximum and falls off for larger k. As 
in the previous examples, the leading order of the uni- 
form approximation has an amplitude error of roughly 
10% with respect to the exact numerical results. The 
(second-order) improved leading order uniform approx- 
imation, however, lies almost on top of the numerical 
results. Remarkably, although the shape deviates from 
a simple power-law behavior quite significantly, the im- 
provement strategy is still effective. 

The relative errors are shown in Figs. 113b and 113b . 
The error behavior divides into two regimes, to the left 
and the right of k ~ 0.05 hMpc^ 1 . The behavior to 
the right is that of a 2 -model [Cf. Eqn. while 
the behavior to the left is that of a polynomial poten- 
tial with linear, quadratic, and quartic terms (|70|l . Note 
that the error estimate from Eqn. (|47|l changes sharply 



across this divide, by more than an order of magnitude, 
from - 0.1% to ~ 0.005%. To the left, this error esti- 
mate shows that there is no point in attempting a cor- 
rection beyond second-order using Eqn. I|4fj[l , consistent 
with the results shown for second and fourth-order cor- 
rected spectra. To the right, the smallness of the error 
estimate is consistent with the improved quality of the 
fourth-order results. Results for the tensor spectrum are 
qualitatively similar. As expected, the uniform approx- 
imation improves on the slow-roll result to the left of 
k ~ 0.05 /iMpc -1 , since the slow-roll assumptions are 
violated in this region. 

The spectral indices are displayed in Figs. 114b and!14b. 
The potential z" / z (see Fig. Illf) leads to a blue scalar 
spectrum for smaller momenta and a red scalar spectrum 
for larger momenta (Fig. 114b ). For the spectral index 
the uniform approximation in leading order is remarkably 
close to the exact numerical result, more or less indepen- 
dent of k. It only deviates slightly at k ~ 0.04 ft, Mpc -1 
for the scalar spectral index. In Fig. ll5l the relative errors 
of the uniform approximation and the slow-roll approxi- 
mation are displayed. Away from the transition fc-value, 
the relative error is smaller than ~ 0.2% for the scalar 
spectral index and smaller than ~ 0.5% for the tensor 
spectral index. The slow-roll approximation, by compar- 
ison, deviates by ~ 2% from the exact numerical results. 

In all the computations so far, the background equa- 
tions were not approximated to obtain H(t) and <f>(t) and 
its derivatives, but were solved numerically. Based on a 
Taylor expansion in t he p otential V{4>) and its deriva- 
tives, Liddle and Lyth |50| introduced an approximation 
for the background equations. This approach (formal- 
ized and expanded to higher orders by Liddle, Parsons, 
and Barrow |43|) leads to a simplification of the slow-roll 
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FIG. 15: Relative error of the (a) scalar and (b) tensor spectral index (see Figs. 114^ , and !14b ) for the C 2 -potential in Eqns. I69H 
and J7DJ. The green band is the error estimate for the uniform approximation. In addition to the uniform and the slow-roll 
approximation we also show the slow-roll approximation including higher order derivatives of V((f>) [see Eqns. j^D an d Q72JI]. 
Light brown dotted line: first order slow-roll, dark brown dashed-dotted line: second order slow-roll. 
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FIG. 16: Ratio R(k) of tensor to scalar perturbations for 
the C 2 -potential in Eqns. 1691 and 1701 : the green band is 
the error estimate for the uniform approximation to leading 
order; the relative difference between all three approximations 
is below 2% (see Fig. 1171 . making it difficult to distinguish the 
curves. 



FIG. 17: Relative error for the ratio R(k) (ratio of tensor to 
scalar perturbations for the C 2 -potential in Eqns. 1691 and 
fTTijl : the green band is the error estimate for the uniform 
approximation to leading order, the darker green band is the 
estimate for the improved leading order. 



parameters in the following form: 

1 fV'\ 2 _ 1 fV^\ 4 1 V' 2 V" 



Si 



2\V 



i fv'\ 2 v" 2 4 



2 V V 



V 3\V 
4 V' 2 V' 



(71) 



(72) 



3 V V J 3 V 3 ' 
Here we have followed the conventions of Stewart and 



Lyth |39|. Being an asymptotic expansion, this approxi- 
mation can sometimes lead to extra errors. For the last 
model, we have calculated the spectral indices with this 
additional approximation and the degradation in relative 
error is shown in Fig. 1151 

The ratio R(k) of tensor to scalar perturbations is de- 
picted in Fig. while the corresponding relative errors 
for the different approximations are shown in Fig. 1171 
The uniform and the slow-roll approximation are both 
quite close (~ 1% error) to the exact numerical result, 
even though the variations in R(k) are not small. Follow- 
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ing Section IIII Dl the accuracy of the leading-order uni- 
form approximation for the ratio R(k) can be improved 
using Eqn. (|46|1 and the corresponding equation for ten- 
sor perturbations. We have not displayed this improved 
ratio in Fig. 1161 since it would be almost indistinguishable 
from the exact numerical result. However, Fig. 1171 shows 
that the relative error of the second-order improved lead- 
ing order of the uniform approximation is smaller than 
~ 0.3% over the whole /c-range, consistent with the error 
estimate (|T7J) . 



D. Consistency relations 

Single field inflation is characterized by degeneracies in 
the observable parameters, such as the tensor to scalar 
ratio R(k), the scalar spectral index ns, the tensor spec- 
tral index ny, and their respective running. Not all the 
parameters are independent; they are connected by so- 
called "consistency relations." For power-law inflation, 
e.g., the tensor spectral index is related to the scalar 
spectral index via the special result tit = ns — 1- In the 
context of slow-roll inflation, consistency relations be- 
tween, e.g., R(k) and nx{k) have been derived (see e.g., 
Ref. |5fl | for an early derivation, note, however, that Lid- 
die and Lyth define R{k) as the ratio of the quadrupole 
moments. Since this definition introduces further uncer- 
tainties through the transfer functions here we use the 
amplitudes of the power spectra themselves). A more 
general expectation beyond slow-roll is that R(k) is some 
function of ns{k) and nr(k). 

Following Ref. we introduce rescaled amplitudes 
via 



0.6 



A s (k) 
A T (k) 



2 p y;i 
5 



P%'(k), 



1 

10 



(73) 
(74) 



leading, in the first order slow-roll expansion, to a con- 
sistency relation between R(k) and nx(&) of the form: 



1%T 



J2 
, /i T 



Or, adapted to our notation, 

P h (k) 8P T {k) 



R(k) 



Pn(k) P s (k) 



-8n T (k). 



(75) 



(76) 



The equality in this relation is only valid when the first- 
order slow-roll quantities are used, even for power-law 
inflation. Note also that the relation is fc-dependent in 
general. 

Inspired by the idea of consistency relations between 
measurable quantities from the microwave background, 
Dodelson et al. [U classified inflationary models depend- 
ing on the relative magnitudes of the spectral indices ns 
and tit, and the ratio of the amplitudes of the power 
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FIG. 18: Tensor to scalar ratio R(k) against the scalar spec- 
tral index ns(k) for the exact numerical results of the inflation 
models investigated here; via the k dependence (running of 
the spectral index and the tensor to scalar ratio) each model 
produces a line - the solid red line denotes the <^ 2 -model, the 
dashed green line the </> 4 -model and the dashed-dotted blue 
line the glued C 2 -potential; the black crosses denote the cen- 
tral WMAP pivot k = 0.05 Mpc -1 . The two dotted lines 
delimit small field, large field and hybrid inflation models ac- 
cording to the classification in Ref. |51f . In the insert we show 
a zoom into the region where the C 2 -potential goes over into 
the 2 -potential. It is interesting to note that at this point 
the curve bends over sharply. 



spectra R(k). In terms of slow-roll parameters, the clas- 
sification is based on the relative magnitudes of e and S± . 
Using Eqn. (|76|l . relations between slow-roll parameters 
can be translated to inequalities between R(k) and ns- 
The three model types are then characterized by: 



(i) Small field models: 

-Si < -e or R < -§(n<? - 1), 

(ii) Large field models: 

-e < 5i < e or -|(n s - 1) < R < -t 

(iii) Hybrid models: 

e < -tfi or -8(n s -1)<R. 



(n s -l), 



In Fig.^Jwe show R(k) as a function of ns(k) for the 
three models investigated in the previous sections. The k- 
dependence, i.e., the running of the quantities, produces 
lines, rather than single points in the R — ns plane. For 
the single-field chaotic models these lines are straight. 
This is not the case, however, for the C 2 -potential where 
the fc-dependence is not a straight line, but a curve. For 
all three models we have marked the location of the cen- 
tral pivot fc* = 0.05 Mpc -1 with a cross. In all cases the 
total number of e-folds is fixed to be N ~ 60 counted 
from the pivot scale. 

We can also test the consistency relation R ~ —871^ of 
slow-roll inflation. In Fig. ^3 R(k) is plotted as a func- 
tion of nr(k). The slow- roll dominated chaotic models 



21 



1 1 1 1 1 1 — 2 1 

"X*— k=10 5 /Mpc <t> 4 

(S 




glued V( 0) 

k=0.05/Mpc + 

-8n T (k) 







-0.040 -0.035 -0.030 -0.025 -0.020 -0.015 -0.010 -0.005 0.000 
n T (k) 

FIG. 19: Tensor to scalar ratio R(k) against the tensor spec- 
tral index nrik); specifications as in Fig. 1181 except that the 
dotted line here denotes R = — 8nr- By fitting the <f> 2 - and 
the 4 -model results, they are found to lie on a straight line 
R = —(7.757 ± 0.004)riT. Although the numerical prefactor 
is different from the slow-roll estimate of —8, the degeneracy 
is still obvious. As in the previous figure, the insert shows 
a zoom into the region where the curve of the C 2 -potential 
bends over to join the line from the (^-potential. 



are indeed very close to the slow-roll consistency rela- 
tion R ~ — 8riT (the dotted line in Fig. 1191) . but not the 
C 2 -model. A detailed understanding of the consistency 
relations and their observational value requires further 
study. 

Regarding observations, the allowed range in the R- 
ns plane is already constrained by WMAP, the SDSS 
Galaxy clustering information and the SDSS Lya for- 
est measurements (see e.g., Ref. H§[). An observation 
of parametric degeneracy by future CMBR experiments, 
e.g., PLANCK 52], would provide strong support for in- 
flation as the source of primordial fluctuations. 



VII. CONCLUSIONS 

In this paper we have presented different approaches 
to the computation of primordial power spectra and the 
corresponding spectral indices and their running, with a 
view to understanding and controlling the various sources 
of error in the calculations. In addition, the ratio of the 
tensor to scalar power spectra and the related consistency 
relation were investigated. 

We have implemented an efficient and accurate method 
for exact mode-by-mode integration utilizing results from 
the uniform approximation to set up the initial condi- 
tions. The uniform approximation, introduced for calcu- 
lating inflationary perturbations in 0, 0] , was numeri- 
cally investigated in detail. We showed that the leading 
order results in this approximation for the power spec- 
trum can be easily improved for well-behaved i>s and vt 
using previously obtained results for the case of constant 



v (the results for the spectral indices n$ and tit are al- 
ready excellent and do not need to be improved). Our 
numerical and semi-analytic results for the power spec- 
tra and spectral indices agree within 0.1%. Thus, the 
primordial power spectra can be determined at the same 
level of accuracy as the transfer functions. 

A key feature of the uniform approximation is the ex- 
istence of an error control theory, which is missing in the 
slow-roll approximation. At leading-order in the uniform 
approximation, we showed how to implement a useful 
approximate error bound for the power spectrum. In 
addition, we provided an error estimate for the power 
spectra obtained from the improved first order uniform 
approximation and the spectral indices. These error es- 
timates are based on our knowledge of the ultra-local 
corrections from the second order uniform approxima- 
tion results. Thus, these estimates are much tighter than 
the general leading order error bounds from the uniform 
approximation. 

We have analyzed in detail three different classes of 
inflationary models: power-law models, chaotic models, 
and a model with a C 2 -potential. We used the power- 
law model, for which exact analytical results exist, to 
demonstrate the accuracy of our numerical implementa- 
tions of the exact mode-by-mode integration, the uniform 
approximation, and the slow-roll approximation. The de- 
viation of the analytic and numerical results was in gen- 
eral very small, around 1 part in 10 6 . Two scale- free 
chaotic inflationary models were chosen as representa- 
tives for common slow-roll models: For these models the 
slow- roll approximation was excellent. As a final exam- 
ple, we constructed a C 2 -potential which had two dy- 
namical phases patched together in a relatively smooth 
manner. In the first phase, where v changes more rapidly 
than in the second phase, the uniform approximation was 
much more accurate than the slow-roll approximation, 
while in the second phase both approximations produced 
very good results with small errors. 

Up-coming high-precision CMBR measurements will 
provide data to constrain the zoo of inflation models. An 
accurate and fast code for calculating primordial power 
spectra, spectral indices, and their running will be crucial 
to this analysis. We are developing an interface to con- 
nect our code to Boltzmann solvers |5i| in order to gen- 
erate the C;'s directly. In addition, we will use the code 
to test the robustness of the information obtained on the 
inflationary equation of state from the measured power 
spectrum by our recently introduced non-parametric re- 
construction program |j| ■ These two complementary ap- 
proaches provide new precision probes of the first mo- 
ments of our Universe. 
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APPENDIX A: NUMERICAL 
IMPLEMENTATION DETAILS 

1. Comment on the Momentum Discretization 

The actual momentum discretization chosen for nu- 
merical work is arbitrary but it is a good idea to adjust 
the chosen values of the momenta so that the relation 



9 2 {f)) = k 2 f => <?(M) = o, 



(Al) 



defining the momentum-dependent turning points fj — 
fj(k) is satisfied exactly, even though the conformal time 
is known only at discrete points. This can be done by 
locking the momentum discretization to the time dis- 
cretization, i.e., by guaranteeing that if the time dis- 
cretization is given, k discretization points are chosen 
only if they satisfy Eqn. (|A1|) . Of course a predefined 
momentum discretization is unnecessary if we only wish 
to calculate power spectra and spectral indices in the 
uniform approximation; the predefined momentum dis- 
cretization is used only for initializing the mode func- 
tions in the mode-by-mode approach, where we need the 
integrals on the left of the turning point. 



2. Spectral Indices in the Uniform Approximation 



chosen. It is further required that —2fj > Arj, i.e., Ary 
be sufficiently small. As f\{k) — > CP in the limit k — oo, 
this relation also constrains the highest reliable mode for 
a given time discretization in the exact numerical results. 
The first integral gives 



ri+Ar/ 



drj' 



Vgik^f) 



drj 



k 2 



= -^-2rjA v - Aif. (A3) 

In order to avoid calculating the integrals numeri- 
cally up to r\ — > CP, we calculate the remainder of the 
integral from an asymptotic value r] a , where the inte- 
grand is sufficiently small and we can stop the numer- 
ical integration, to 77 ~ CP, assuming that v 2 (r\) ~ 
v 2 (r) a ) + 2v{r\ a )v' (r\^){r\ - r/ a ). Then we have 



-2k 2 lim 



dr]' 



1?*- 



k 2 



-2^J v 2 (r) a ) - 2v{r] a )v' (r/ a )?7 a + 2i/^ 2 (7?a) - k 2 r\i 
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-arcsm- 



'A 

1v{r) a )v '(77 a 



A 



,(A4) 



with 



A = [2v^ a )v\^)\ 2 +±k 2 [v 2 (r ]a )-2 V ^ a )v\^)"q a \. (A5) 

Alternatively, it is possible to convert the evaluation 
of the spectral index into the problem of solving a differ- 
ential equation, rather than evaluating an integral. The 
limit ki] — > CP in Eqn. I|48(l is interchangeable with a con- 
formal time derivative, so that the physical time deriva- 
tive of the spectral index reads 



The integral for the spectral index has a square root 
singularity at 77 — fj and is handled specially in the nu- 
merical routine. We split the integral appearing in the 
spectral index into two parts: 



di 



rj+Ar] 



drf 



drj 



\/g(k,v') Jfj+A V \AKfc,r?')' 

(A2) 



where A77 is a small quantity. Note that At7 is really k 
dependent, because the discretization in 77 is not equidis- 
tant. 

In the first integral we can substitute v 2 (j\) by v 2 {ff), 
i.e., insert the leading order of the local approximation. 
The first integral can then be calculated analytically and 
keeps track of the inverse square root singularity, while 
the second integral has no singularity and can be easily 
calculated numerically. The quantity Arj is given by the 
actual time discretization in physical time t that we have 



ns[k,ri(t)} 



-2k 2 



a(t)y/gs(r],k) 



(A6) 



In order to avoid the square root singularity the integra- 
tion starts at 77 = 77 + A77, so that the "initial" condition 



ns 



(fc, 77 + A77) = 4 - 2k\/ -2f] Arj - Aff 



(A7) 



includes the integral in Eqn. 1A3|) . It is understood that 
the limit kr\ — * CP is taken when calculating ns(k). The 
integration of the differential equation with a high order 
integrator is more precise than a standard trapezoidal 
rule. In contrast, as the discretization in 77 is not equidis- 
tant, higher-order integration schemes would be some- 
what more complicated to implement. However, we have 
verified that a standard trapezoidal integration rule al- 
ready gives sufficiently precise answers. In fact, the local 
approximation for the spectral index at leading order is 
quite close to the numerical nonlocal integral in the cases 
where the derivative expansion is valid. 
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3. Conversion to Physical Units 

For completeness, we explain here how units are han- 
dled in the numerical implementation. As always, it is 
convenient to work in dimensionless units, i.e., by choos- 
ing H = c = 1. In addition, we set the factor 8irG in 
the Fricdmann equation in the numerical code to unity. 
These choices lead to values for the input parameters, 
e.g., initial conditions and coupling constants, of order 
unity. This helps prevent numerical problems arising 
from very large or very small numbers. In order to recon- 
vert the dimensionless units to physical units the Hubble 
parameter H has to be rescaled via 



(A8) 



where H is the dimensionless Hubble parameter used in 
the code and Hq = lOOh kms Mpc . The rescaled 
momentum k in physical units h Mpc -1 is therefore given 
by 



^phys 



100ca(0) 



k kms x Mpc \ 



(A9) 



with c = 2.99792458 x 10 5 kms -1 . The initial expansion 
rate is a(0). Throughout the paper we have dropped 



the suffix "phys" implying that all results are given in 
physical units. 



Next we discuss the normalization of the amplitude of 
the power spectrum. The power spectra for scalar and 
tensor perturbations as defined in Eqns. Ijlll) and 112(1 are 
dimensionless and therefore not sensitive to the units of 
k. Their amplitude is determined fully by the parameters 
chosen in the inflaton potential V (</>). Since parameters 
such as the inflaton mass m 2 in a chaotic m 2 2 -model 
are generally not known, we present the results for the 
power spectra with respect to the WMAP normalization 
where the amplitude of scalar perturbations is given by 
|Ai? 2 | = 2.95" 9 ^ with A = 0.9±0.1 (at K = 0.05/Mpc), 
see Refs. ^H, ^J. Using the fact that for a fixed num- 
ber of e- folds, counted from the horizon crossing of 
the parameters in the monomial potentials simply lead 
to a global normalization factor [24| we can avoid hav- 
ing very small numbers in the numerical calculations and 
just normalize the spectra afterwards. When stating re- 
sults we give in each case the parameters in the potential 
corresponding to the WMAP normalization. 
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